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We study the order parameter symmetry in a hexagonal crystal with co-existing superconductivity 
and ferromagnet ism. An experimental example is provided by carbon-based materials, such as 
graphite-sulfur composites, in which an evidence of such co-existence has been recently discovered. 
The presence of a non-zero magnetization in the normal phase brings about considerable changes in 
the symmetry classification of superconducting states, compared to the non-magnetic case. 
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The recent discoveries of superconductivity co-existing 
with ferromagnetism in such transition element com- 
pounds as ZrZn 2 and UGe 2 || have called for a re- 
vision of our views on the interplay of the two phenom- 
ena. Symmetry considerations can help us to identify the 
possible order parameters, even if the microscopic mech- 
anism of pairing is not firmly established, which is often 
the case for the systems with strong electron correlations. 
The standard group-theoretical classification of the su- 
perconducting states in crystals || (for a review, see Ref. 
j|) is not applicable if the normal state is ferromagnetic. 
In Refs. ||, the symmetry analysis was extended to cover 
the magnetic case, and also a complete classification of 
the superconducting states in a cubic ferromagnet such 
as ZrZn 2 was done. The orthorhombic symmetry, which 
is appropriate for UGe2, was studied in Refs. j6|. 

Very recently, a strong evidence was found of super- 
conductivity apparently co-existing with ferromagnetism 
in yet another class of materials, namely graphite-sulfur 
composites (see also Ref. ||). The superconducting 
and ferromagnetic behavior in highly-oriented pyrolitic 
graphite was reported in Ref. || . The experimental data 
suggest that the superconductivity is associated with the 
graphite planes, whose symmetry is hexagonal. In this 
article, we study the possible pairing symmetries and the 
related nodal structures of the superconducting gap in a 
hexagonal ferromagnet, thus complementing the analy- 
sis in Refs. The details of the electronic spec- 
trum of graphite [ flO| |, the nature of ferromagnetic and 
superconducting correlations, and also the fact that its 
single layer has semimetallic properties with some pe- 
culiar physics related to that do not influence our 
results. We only assume that a standard BCS-Gor'kov 
phenomenology is applicable, i.e. the origin of supercon- 
ductivity is the Cooper pairing of fermionic quasiparticles 
with opposite momenta, so that the order parameter is 
an anomalous average of the pair creation operator. 

The starting point of the group-theoretical analysis 
is the symmetry group Q of the normal state which is 
defined as a group of transformations which leave the 
system Hamiltonian invariant. In non-magnetic super- 
conductors, time reversal symmetry is not broken, and 
Q = S x K x f/(l), where S is the space group of the 



crystal, K is time reversal operation, and U(l) is the 
gauge group. In contrast, in magnetic crystals time re- 
versal symmetry is broken, and Q = Sm X U(l), where 
Sm is the magnetic space group which is a group of sym- 
metry operations leaving both the crystal lattice (the mi- 
croscopic charge density) and the magnetization density 
M invariant jl2|. For example, if there is a crystal point 
group rotation R which transforms M to — M, then the 
combined operation KR will be an element of Sm, be- 
cause the time reversal restores the original M not affect- 
ing the lattice symmetry. We assume that the spin-orbit 
coupling is sufficiently strong, so that the space group 
elements act on the orbital and spin coordinates simulta- 
neously. As shown in Refs. fl^j , the symmetry analysis 
should be modified in the magnetic case, due to the fact 
that the operation KR is anti-unitary. 

The crystal structure of graphite is described by layers 
of honeycomb lattices of carbon atoms, the precise ar- 
rangement of the layers along the z axis being not impor- 
tant for our analysis. We make the usual assumption that 
the subgroup of translations (which are replaced by mag- 
netic translations in our case) does not affect the momen- 
tum dependence of the order parameter. Therefore, the 
space group in the paramagnetic phase can be replaced 
by the hexagonal point group T>eh = D6 x Ci; where 
Cj = {E, 1} (we assume that the space group contains 
the inversion operation /) . The group D 6 is generated by 
the rotations C§ z about the z axis by an angle 7r/3, and 
the rotations Ci x about the x axis by an angle 7r. In the 
ferromagnetic phase, there are three possibilities for the 
easy direction of magnetization: M can be along the z 
axis, the x axis, or the a axis, where a = (v^/2, 1/2,0). 
Here we study only the first possibility; the other cases 
can be treated in a similar fashion. If M is along the z 
axis, then the symmetry of the crystal is described by the 
magnetic group D 6 (C 6 ) x Cj, which is generated by: (i) 
the rotations Cq Zi (ii) the operations KCi x which com- 
bine the rotations Ci x with the time reversal, and (iii) 
the inversion /. The subgroup in parentheses (the unitary 
subgroup) incorporates all symmetry elements which are 
not multiplied by the anti-unitary operation KC2 X , i-c 
D 6 (C 6 ) = C 6 + KC 2x x C 6 . [If M is along the x or 
a axes, then the symmetry is described by the magnetic 
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group D 2 (C 2 ) x C u where D 2 (C 2 ) = C 2 + KC 2z x C 2 .] 
In the presence of spin-orbit coupling, spin is not a 
good quantum number and is replaced by pseudospin. 
The free electron Hamiltonian Hq in the normal state is 

= X)fc[ e +( fc ) c i+ c fc+ + e -( fe ) c fc- c fc-]> where e±(k) are 
the energy spectra for the pseudospin-up and pseudospin- 
down sheets of the Fermi surface. The transformation 
properties of the single-electron states \k,+) and |fe, — ) 
with respect to the magnetic symmetry elements are the 
same as those of the spin eigenstates \k,1) and \k, |), 
which means that Cq z c\ ± C 6z 1 = e T47r / 6 c, + 
{KC 2x ){\cl ± ){KC 2x )- x '= ±i\*d 



-C 2x fc,±' 



C 6i fc,±' 

and 



Ic k ±1 1 = c^ fc± . Here A is an arbitrary c- number 
(note that the combined operation KC 2x is anti-linear). 
The transformation properties of the superconducting 
order parameter can be derived using the mean-field 
approximation for the pairing Hamiltonian: 



' a,(3=± 



(1) 



From the anticommutation relations for c^ a , we see 
that A ++ (fc) and A__(fc) are odd functions of k, but 

A -| (k) = — A |_(— k) does not have a definite parity. 

The order parameter matrix can be cast in a more 
familiar form A(k) = (icra 2 )d{k) + (ia 2 )do(k), where 
d{k) = —d(—k) and do(k) = do(—k) are the pseudospin- 
triplet and the pseudospin-singlet components respec- 
tively . It is convenient to introduce an orthogonal ba- 
sis of unit vectors e%, e 2 , e 3 in the pseudospin space, such 
that e 3 || M and use the following representation of the 
vector order parameter: d{k) — e + cL(fc) + e_d + (fc) + 
e 3 d 3 (k), where e± — (ei ± ie 2 )/\/2 and d + = {d\ + 
id 2 )/V2 = A__ /V2, d_ = (di - id 2 )/V2 = -A++/V2, 
d 3 = (A+_+A_+)/2. Also, d = (A + _-A_+)/2. From 



Eq. (|l|) , we obtain the transformation rules for the order 
parameters under the magnetic group operations: 

C 6z d ± (k) =e ± "/ 3 d ± (C^ z 1 k) 
C 6z d 3 (k) - daiC^k) 
KC 2x d ± (k) = dUC^k) (2) 
KC 2x d 3 (k) = -d* 3 (C^k), 

and 

C 6z d (k) = d (Ce z 1 k) 
KC 2x d (k)=d* Q (C 2x 1 k). (3) 

Because of the inversion symmetry, the triplet and the 
singlet order parameters can be considered separately. 

Triplet order parameter. The superconducting order 
parameter transforms according to one of the irreducible 
representations T of the normal state symmetry group 
Q. In our case, Q contains the anti-unitary operation 
KC 2x , and the standard symmetry analysis should be 
modified: instead of usual representations, one should use 
co-representations of the magnetic point group De(Cg), 



which can be derived from one-dimensional representa- 
tions of the unitary subgroup Q. The results for 
odd co-representations T are listed in Table S. Note that 
the action of the unitary and anti-unitary orbital symme- 
try elements on scalar functions f(k) is defined as follows: 
Rf(k) = /(iT 1 *;), and KRf{k) = f^-R^k). 

The basis functions of the co-representations in Table 
[| vanish on some lines or planes in the momentum space. 
Although the specific form of the basis function is not 
imposed by symmetry [for example, any fr(k) can be 
multiplied by an arbitrary real function of k^, k\ + ky, 
(k y ± ik x ) 6 , etc.], the position of the zeros in fc-space is 
independent of the choice of fr{k). It can be checked 
that the zeros of the basis functions from Table | are all 
required by symmetry, and, on the other hand, there are 
no other symmetry-imposed zeros of the basis functions. 

All co-representations of D6(C@) are one-dimensional, 
therefore the order parameter can be represented as 

d r (fc) = iri-e + f T _(k) + i-q + e-f r+ (k) + r? 3 e 3 /r 3 (fe). 

(4) 

Thus, the order parameter has three components rj+, 
?7_, and r/ 3 , which enter the Ginzburg-Landau (GL) free 
energy and can be r-dependcnt, in general. It is ob- 
vious from Eqs. (0) and (|J) that the orbital symme- 
tries of d+, d_, and d 3 , which are described by the co- 
representations r + , r_, and r 3 respectively, should all 
be different. In Table we list the combinations of 
the orbital co-representations which give rise to the same 
symmetry T of the vector d. For example, the vector 
d(k) = ir}-e+f Bu (k) + ir} + e-fx Elu (k) + ij 3 e 3 f2 E2u (k) 
transforms according to the co-representation 2 E 2u . 

The transformations of the pseudospin vector d under 
the magnetic symmetry elements can be interpreted as 
operations acting on r\ = ?y_, rj 3 ) T ■ Namely, for a 
co-representation T, Ce z t](r) — Xr[Ce z ]r](r), where xr 
is the character of Cq z from Table |. Also, because of our 
choice of the overall phases of the basis functions (see the 
caption to Table |jj) and the presence of the factors i on 
the right-hand side of Eqs. (Q), KC 2x r](r) = rj* (r). Note 
that there are no symmetry operations that transform, 
say, r)+ to ??_, etc, which has profound consequences for 
the GL theory of ferromagnetic superconductors. Indeed, 
the GL free energy contains all combinations of the or- 
der parameter components and the spatial gradients that 
are invariant under the transformations from Q . The fact 
that the only change that occurs to r\ under the magnetic 
symmetry operations is the multiplication by a phase fac- 
tor which is the same for all three components means 
that there are more invariant terms allowed in the GL 
functional for a ferromagnetic superconductor than for 
its non- ferromagnetic counterpart. The general form of 
the uniform terms in the free energy density is 

F= AijiT^rjj + Bij,ki'n*rf j r) k 'n u (5) 

ij=±,3 ijkl=±,3 

where Ay is a real symmetric matrix, and the matrix 
B is symmetric with respect to i <-> j and k <-» I, and 
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TABLE I: The character table and the examples of the 
basis functions for the odd irreducible co-representations of 
the magnetic point group De(Ce). The overall phases of 
the basis functions are chosen so that KC2 X fr(k) = fr(k). 
uj = exp(27ri/3), and Ai,2 are arbitrary real constants. 



r 


E 




fr(k) 


A u 


1 


1 


k z 


B u 


1 


-1 


\i(k y + ikx) 3, + X 2 (k y — ikx) 3. 


1 Eiu 


1 


— LU* 


ky ~\~ ikx 


2 Eiu 


1 


—LU 


ky ikx 


x E 2 u 


1 


LU* 


kz{ky ikx^) 


2 E 2u 


1 


LU 


kz^ky -\- ikx^) 



TABLE II: The sets of orbital co-representations correspond- 
ing to the same symmetry of the triplet order parameter dr(k) 
[see Eq. (§)]. 



r 


(/r + (fc),/r_(fc),/r 3 (fe)) 


Au 


( 2 Eiu, 1 Eiu, Au) 


B u 


( 2 E 2u , 1 E 2 u, B u ) 


1 Eiu 


(A u , 2 E 2 u, 1 Eiu) 


2 Eiu 


( 1 E 2u , A u , 2 Eiu) 


1 E 2u 


(B u , 2 Eiu, 1 E 2u ) 


2 E 2u 


( 1 Eiu, B u , 2 E 2u ) 



satisfies the following condition: Bij^i — Bki^j. The 
critical temperature T c is defined as the temperature at 
which one of the eigenvalues of A changes sign. At T > 
T c , A is positive definite, and r)+ = r/_ = rj3 = 0. Below 
T c , all three components of rj are non-zero, in general. In 
addition to Eq. (Q) , the GL functional contains a variety 
of the gradient terms. All this can lead to a rich phase 
diagram, which we shall not discuss here. Let us just 
note Eq. (|J) is formally equivalent to a model of a three- 
band superconductor with three scalar order parameters 
of the same symmetry. 

An important consequence of the above results is that 
the gap in the spectrum of Bogoliubov quasiparticles goes 
to zero at some points or lines at the Fermi surface. The 
excitation spectrum can be obtained by diagonalizing the 
Hamiltonian H = Hq + Hmf- The quasiparticle energy 
E(k) vanishes at some k if the following condition is sat- 
isfied: 



2e A 



+ dl\ 2 = 0. 



(6) 



Thus, the condition for the gap zeros on the pseudospin- 
up sheet of the Fermi surface [at e+(k) = 0], is that 
cL(fc) = d 3 (k) = 0, i.e. / r _(fe) = /r 3 0) = 0. For the 
gap zeros on the pseudospin-down sheet [at e_(fc) = 0], 
we must have d+(k) — dj(k) — 0, i.e. /r + (fc) = fr 3 (k) = 
0. Using Tables | and |J, we see that if the symmetry 
of d corresponds to the co- representations B u , 1 £ , 2u, or 



2 £ , 2„, then all three orbital basis functions have zeros on 
the line k x = k y = 0, so that the energy gap vanishes 
at the poles of both sheets of the Fermi surface. On the 
contrary, for the E\ u and 2 E\ U symmetries, the gap goes 
to zero at the poles of one of the sheets, while the other 
sheet remains fully gapped. For the A u symmetry, there 
are no symmetry-imposed gap zeros. 

In the discussion above, we implicitly assumed that 
all three components of d have comparable magnitude. 
More realistic scenario is that the conditions for the 
appearance of superconductivity are more favorable on 
one of the sheets of the Fermi surface, so that only 
one component of the order parameter, say d- on the 
pseudospin-up sheet, is dominant. Then, the d + compo- 
nent is induced by the inter-band interactions of the form 
c fe+ c -fc + c fc'- c -fc'.-i which are expected to be small if 
the spin-orbit coupling is weak (these interactions vanish 
at zero spin-orbit coupling because of the spin conserva- 
tion). Also, if the exchange band splitting greatly exceeds 
the superconducting energy scale, which is of the order of 
T c , then the inter-band interactions <J k+ c}_ k _Ck'—C—k',+ 
responsible for the ^-component, are negligibly small. 
While this is the case in such materials as ZrZn2 and 
UGe2, it is not clear whether it is also true for the 
graphite-based ferromagnets. If the <i3-component can 
indeed be neglected, then the conditions for the presence 
of the gap zeros become less restrictive, so that, as seen 
from Tables | and |l| the energy gap always has either line 
and/or point nodes on both sheets of the Fermi surface. 

Singlet order parameter. Similar to the ^-component 
of the triplet pairing, the singlet pairing can only be real- 
ized if the exchange band splitting is less than the super- 
conducting T c (Chandrasekhar-Clogston limit) fl4| . This 
limit can be slightly exceeded if to consider the possibil- 
ity of the Cooper pairs having a non-zero momentum 
(Larkin-Ovchinnikov-Fulde-Ferrell state) |l5fl . The sym- 
metry analysis can be done similarly to the triplet case. 
The only difference is that there is only one order param- 
eter component, and the symmetry is described by even 
co-representations of Dg(C6). We have 



do,r(fe) = ^/r(fe), 



(7) 



where -0 is a quantity which enters the GL functional, and 
/r(fc) is the basis function for the co-representation T, see 
Table pl[ As mentioned above, fr(k) can be multiplied 
by an arbitrary real function of k\, k% + ky, (k y ± ik x ) e , 
etc. Under the action of the magnetic symmetry ele- 
ments, C% z xp(r) — xr[C6z]ip(r), where xr is the charac- 
ter of C Gz from Table HI, and also KC 2x 4'( r ) = ^{i")- 



The condition for the gap in the excitation energy E(k) 
to vanish at some k is simply 



e + e_ + \d \ 



0. 



(8) 



therefore the gap nodes appear simultaneously on both 
sheets of the Fermi surface where the basis fun ction /r(fc) 
has symmetry-imposed zeros. From Table III, we see that 
for all order parameter symmetries, except from A„, the 
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TABLE III: The character table and the examples of the basis 
functions for the even irreducible co-representations of the 
magnetic point group D6(Ce). For all F, KC2xfr(k) = fr(k), 
and Ai,2 are arbitrary real constants. 



r 


E 




k(k) 




1 


1 


1 


B a 


1 


-1 


k z [\i{k y + ik x f + \ 2 {k y - ik x f] 


Elg 


1 


— LU* 


k z {ky -\~ ik x ) 


2 Elg 


1 


— LU 


kzi^ky ik x ^) 


E 2 g 


1 


LU* 




2 E 2 g 


1 


LU 


(ky + ik x ^ 



gap goes to zero either on the equators or at the poles of 
both Fermi surfaces. 

The expressions for the basis functions of the co- 
representations of Dg(Ce) given in Tables [l| and [II 



are 

applicable if the Fermi surface is centered around the 
r point of the first Brillouin zone. It is not the case 
in graphite, where the Fermi surface consists of small 
"sausage" -like pockets along the six vertical edges of the 
hexagonal Brillouin zone, i.e. along the lines k = Cg z k s , 
where n — 0, ...,5, k s — 2Ki/3 — K-z/Z + k z e z , and 
Ki and K2 are the reciprocal lattice vectors fnJ. It 
can be easily proved that for all odd and even E co- 
representations, in addition to the zeros at k x — ky — 0, 
the basis functions also vanish at the vertical edges of 



the Brillouin zone, because k s and C% z k s are equivalent 
points. Similarly, if a basis function vanishes on the plane 
k z = 0, then it should also vanish at the horizontal sur- 
faces of the Brillouin zone, i.e. at k z = ±7r/co, where cq 
is the lattice constant of graphite along the z axis. 

The presence of the gap nodes would manifest itself 
in power-law temperature dependences of the thermody- 
namic and transport properties [Q. For example, the 
electronic specific heat at low temperatures should be 
C(T)/T ~ T 2 for the first-order point nodes, C(T)/T ~ 
T for the line nodes or the second-order point nodes, and 
C(T)/T ~ T 2 / 3 for the third-order point nodes. 

To summarize, we have derived a complete symmetry 
classification of the superconducting states in a hexag- 
onal ferromagnet, considering both the triplet and the 
singlet pairing channels. We have discussed the nodal 
structure of the gap in the excitation spectrum, and also 
the modifications to the Ginzburg-Landau theory for fer- 
romagnetic superconductors. So far, the only experimen- 
tal system to which our results might be applicable is the 
graphite-sulfur composite studied in Ref. [0. 
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vic for useful comments and interest to this work, and 
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the experimental situation. The financial support from 
Brock University is gratefully acknowledged. 



[1] 
[2] 

[3] 

[4] 
[5] 

[6] 
[7] 



C. Pfleiderer, M. Uhlarz, S. M. Hayden, R. Vollmer, 

H. v. Lohneysen, N. R. Bernhoeft, and G. G. Lonzarich, 
Nature 412, 58 (2001). 

S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, 
R. K. W. Haselwimmer, M. J. Steiner, E. Pugh, 

I. R. Walker, S. R. Julian, P. Monthoux, G. G. Lonzarich, 
A. Huxley, I. Sheikin, D. Braithwaite, and J. Flouquet, 
Nature 406, 587 (2000). 

G. E. Volovik and L. P. Gor'kov, Sov. Phys. - JETP 
61, 843 (1985) [Zh. Eksp. Teor. Fiz. 88, 1412 (1985)]; 
K. Ueda and T. M. Rice, Phys. Rev. B 31, 7114 (1985); 
E. I. Blount, Phys. Rev. B 32, 2935 (1985). 
V. P. Mineev and K. V. Samokhin, Introduction to Un- 
conventional Superconductivity (Gordon and Breach, The 
Netherlands, 1999). 

M. B. Walker and K. V. Samokhin, Phys. Rev. Lett. 
88, 207001 (2002); K. V. Samokhin and M. B. Walker, 
Phys. Rev. B 66, 024512 (2002); K. V . Samokhin and 



M. B. Walker, |arXiv:cond-mat/0206487[ 
V. P. Mineev. |arXiv:cond-mat/0204263i I. A. Fomin, 
arXiv:cond-mat/020715"2| (see also I. A. Fomin, JETP 



Lett. 74, 111 (2001) [Pis'ma Zh. Eksp. Teor. Fiz. 74, 116 
(2001)], for an earlier attempt of the symmetry analysis 
for UGe 2 ). 

S. Moehlecke P.-C. Ho, and M. B. Maple, |arXiv:cond" 
mat/0204006[ 



[8] R. Ricardo da Silva, J. H. S. Torres, and Y. Kopelevich, 

Phys. Rev. Lett. 87, 147001 (2001). 
[9] Y. Kopelevich, P. Esquinazi, J. H. S. Torres, and 
S. Moehlecke, J. Low Temp. Phys. 119, 691 (2000). 

[10] M. S. Dresselhaus and J. G. Mavroides, IBM J. Res. De- 
velop. 8, 262 (1964); S. J. Williamson, S. Foner, and 
M. S. Dresselhaus, Phys. Rev. 140, A1429 (1965). 

[11] J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Phys. 
Rev. B 63, 134421 (2001); D. V. Khveshchenko, Phys. 
Rev. Lett. 87, 246802 (2001); G. Baskaran and S. A. Ja- 
fari, Phys. Rev. Lett. 89, 016402 (2002). 

[12] L. D. Landau and E. M. Lifshitz, Electrodynamics of 
Continuous Media (Pergamon Press, London, 1960). 

[13] E. Wigner, Group Theory and Its Applications to the 
Quantum Mechanics of Atomic Spectra (Academic Press, 
NY, 1959); C. J. Bradley and A. P. Cracknell, The Math- 
ematical Theory of Symmetry in Solids (Clarendon Press, 
Oxford, 1972); A. P. Cracknell, Progr. Theor. Phys. 35, 
196 (1966). 

[14] B. S. Chandrasekhar, Appl. Phys. Lett. 1, 7 (1962); 

A. M. Clogston, Phys. Rev. Lett. 9, 266 (1962). 
[15] A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. - JETP 

20, 762 (1965) [Zh. Eksp. Teor. Fiz. 47, 1136 (1964)]; 

P. Fulde and R. A. Ferrell, Phys. Rev. 135, 550 (1964). 



